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Abstract. Merging binary systems consisting of two collapsed objects are among the most promising sources of 
high frequency gravitational wave, GW, signals for ground based interferometers. Double neutron star or black 
hole/neutron star mergers are also believed to give rise to short hard bursts, SHBs, a subclass of gamma ray 
bursts. SHBs might thus provide a powerful way to infer the merger rate of two-collapsed object binaries. Under 
the hypothesis that most SHBs originate from double neutron star or black hole/neutron star mergers, we outline 
here a method to estimate the incidence of merging events from dynamically formed binaries in globular clusters 
and infer the corresponding GW event rate that can be detected with Advanced LIGO/ Virgo. In particular a 
sizeable fraction of detectable GW events is expected to be coincident with SHBs. The beaming and redshift 
distribution of SHBs are reassessed and their luminosity function constrained by using the results from recent 
SHBs observations. We confirm that a substantial fraction of SHBs goes off at low redshifts, where the merging 
of systems formed in globular clusters through dynamical interactions is expected. 
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1. Introduction 

' Merging binary systems containing two collapsed objects, 
. i.e. a, double neutron star (DNS), a stellar-mass black hole 
' and neutron star (BH-NS) or two stellar-mass black hole 
(BH-BH), are powerful gravitational wave, GW, sources 
' with frequencies from hundreds to over a thousand Hz. 
. They are considered among the most promising GW 
' sources for ground-based interferometers, of the current 
and future generation, such as LIGO, Virgo and their ad- 
vanced versions. BH-NS and, especially, BH-BH mergers 
t emit more powerful and lower frequency GWs than DNS 
mergers, where the sensitivity of LIGO and Virgo detec- 
tors is highest: therefore they can be detected up to larger 
distances. The horizon of first generation LIGO and Virgo 
for DNSs, BH-NS and BH-BH mergers is - 20, 43 and 
100 Mpc, respectively, while Advanced LIGO/ Virgo class 
interferometers should detect them out to a distance of 
~ 300, 650 and 1600 Mpc (for a review see Cutler & 
Thorne 2002). 

The rate of detectable merging events has been es- 
timated based on the observed galactic population of 
DNS binaries containing a radio pulsar (Phinney 1991; 
Narayan, Piran & Shemi 1991; Kalogera et al. 2001; 
Burgay et al. 2003). The best estimate of the DNS merger 
rate in the Galaxy is presently ~ 80tgo°/Myr, convert- 
ing to SOOlgoQ^ Gpc~'^yr~^ for a galaxy number density 



of 10~^Mpc"^ (Kalogera et al. 2004). Population syn- 
thesis studies of binary systems give results consistent 
with the above rate (Perna & Belczynski 2002; Belczynski, 
Kalogera & Buhk 2002; Belczynski et al. 2007). GW sig- 
nals from DNS mergers are expected at a rate of one in 
~ 10-150 years with Virgo and LIGO and one every 1- 
15 days with Advanced LIGO/Virgo class interferometers. 
The BH-NS and BH-BH merger rates in the Galaxy are 
highly uncertain. Belczynski et al. (2007) estimate ~ 1% 
and ^ 0.1% of the DNS merger rate, respectively, implying 
that BH-NS and BH-BH mergers contribute marginally to 
the GW event rate, despite the larger distance up to which 
they can be detected. 

DNSs and NS-BH merging events provide also one of 
the leading models for Short Hard Gamma Ray Bursts, 
SHBs, bursts of 7-rays that last < 2 s and take place at 
cosmological distances (Goodman 1986; Paczyhski 1986; 
Eichler et al. 1989; Narayan, Paczynski & Piran 1992). As 
with the coUapsar scenario for long Gamma Ray Bursts, 
GRBs (Woosley & Mac Fadyen 1999), such models envis- 
age the formation of a black hole surrounded by a torus of 
matter at nuclear density that is rapidly accreted and pro- 
vides the primary source of energy for the burst. For this 
reason it is believed that BH-BH mergers are not among 
the progenitors of SHBs. SHBs comprise about 1/4 and 
1/10 of the CGRO/BATSE and Swift/B AT GRB sam- 
ples. The X-ray and optical afterglows of SHBs have been 
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observed only recently, leading to the first identifications 
and redshift determinations of SHB host galaxies (see e.g. 
Berger et al. 2007a and references therein). Out of the 
~ 30 SHBs in the current Swift and HETE 2 sample about 
half displayed an optical afterglow, 11 have redshift, and 
a few show evidence for beaming in their optical afterglow 
light curve (Fox et al. 2005; Soderberg et al. 2006). The 
association of some SHBs with galaxies characterized by a 
very low star formation rate (SFR, of order <O.1M0 yr~^) 
indicates a long time-delay between the progenitors' col- 
lapse and time at which these SHBs go off. 

It has long been known that DNSs and BH-NS systems 
(hereafter NS-NS/BH systems), can form from massive 
binaries surviving two gravitational collapses (we term 
these "primordial" NS-NS/BH systems). Population syn- 
thesis calculations show that primordial NS-NS/BH sys- 
tems merge after a relatively short time (r w 2 Gyr on 
average), implying that the redshift-distribution of SHBs 
should follow pretty closely the SFR history of the uni- 
verse (e.g. Bloom, Sigurdsson & Pols 1999; Belczynski 
et al. 2002, 2006). Alternatively NS-NS/BH systems can 
form through dynamical interactions in the cores of glob- 
ular clusters, GCs. The dominant process comprises two 
stages: first an isolated NS (or BH) in the cluster captures 
a non-degenerate star during a close encounter, thus form- 
ing a binary containing a collapsed object. This binary 
then undergoes an exchange interaction with another col- 
lapsed object, leading to the ejection of the non-degenerate 
star and formation of a NS-NS/BH binary (Verbunt & 
Hut 1987; Efremov 2000). Grindlay et al. (2006) estimate 
that ~ 10 — 30% of all DNS mergers may stem from 
such dynamically formed systems. Dynamical interactions 
take place efficiently at high stellar densities, which, in 
turn, are highest in the core of globular cluster that have 
undergone core collapse. Therefore, the delay between 
the formation of the collapsed objects and merging of 
dynamically-formed NS-NS/BH systems is dominated by 
the time until core-collapse, which is typically comparable 
to the Hubble time. In their simulations Sadowski et al. 
(2008) found that merging BH-BH binaries form rather 
efficiently in dense clusters, but failed to find any merging 
of NS-NS/BH binaries. 

If NS-NS/BH mergers are among the progenitors of 
SHBs, the SHB rate can be used as an alternative way 
of constraining the DNSs merger rate (Guetta and Piran 
2005, 2006 (GP05, GP06), Nakar et al. 2006). Hopman 
et al (2006) (H06) showed that dynamically formed DNSs 
give rise to more numerous low-z bursts than would be 
expected if SHB closely followed the SFR, as in the case 
of primordial DNSs. The z-distribution of dynamically 
formed DNSs provides a better match to observed red- 
shifts of SHBs. This result was recently confirmed by 
Salvaterra et al. (2008), who also conclude that SHBs may 
well originate from both classes of DNSs. GP05, GP06 and 
H06 fitted the peak flux distribution of SHBs detected by 
BATSE to derive their formation rate and luminosity func- 
tion (LF), in both the dynamically- formed and primordial 
DNS scenarios, while Nakar et al. (2006) and GP06 stud- 



ied the prospects for detecting GW signals from SHBs 
from primordial DNSs. 

In this paper, based on the results from recent SHB 
observations, we first reassess the beaming, redshift distri- 
bution and LF of SHBs. Under the hypothesis NS-NS/BH 
mergers give rise to most SHBs, we then present a method 
to estimate the incidence of merging events from dynam- 
ically formed binaries in globular clusters and infer the 
corresponding GW event rate that can be detected with 
Advanced LIGO/Virgo class interferometers. 

2. Luminosity function and rate evolution of SHBs 

For primordial NS-NS/BH systems, the intrinsic SHB rate 
is given by the convolution of the formation rate of NS- 
NS/BH binaries (which is assumed to follow with negli- 
gible delay the SFR) with the distribution of the merg- 
ing time delays, Ppr(T), where r — tqw is the time over 
which GW losses bring a binary to its pre-merging stage. 
We adopt Ppi-ir) ~ l/r in agreement with the r- values of 
DNS radio pulsar binaries and population synthesis cal- 
culations (Champion et al. 2004; Belczynski et al. 2007). 
The SHB rate from primordial DNSs is given by 

RsHb{z)'^ / dTRsF2{t-T)Pp,{T), (1) 

JQ 

where Rsf2 refers to the "SF2" model for the SFR in 
Porciani & Madau (2001). 

For dynamically formed NS-NS/BH binary, r = rcc + 
'''GWj where the delay time rcc ^ tqw represents the 
elapsed time between the birth of NSs and BHs in GCs 
and the dynamical formation of NS-NS/BH systems fol- 
lowing core collapse. According to H06, the distribution of 
dynamically formed NS-NS/BH binaries Pdyn(''') (which 
replaces Ppr(T) in Eq.l) increases for increasing time de- 
lays. We assume that the formation rate of GCs is pro- 
portional to the total SFR0. We derive an average delay 
of f w 2 and ~ 6 Gyr, for primordial and dynamically 
formed NS-NS/BH system, respectively (for the latter we 
considered only values of r shorter than the Hubble time). 

2.1. The SHB luminosity function 

In this section we summarize the method adopted by 
GP06 and H06 to estimate the SHBs' LF and we derive the 
local SHB rates independently for primordial and dynam- 
ically formed NS-NS/BH binaries. We use here the same 
sample of 194 SHBs detected by BATSE as in Paciesas 
et al. (1999). GP06 and H06 modelled the "isotropic- 
equivalent" LF (i.e. uncorrected for beaming), with a bro- 
ken power law extending from L* / Ai to i*A2 and with 
break at L* (with Ai = A2 = 100), i.e. 

ci>„(L).iogx^co.iog L{ ^y^i^j^,': m 

^ We note that Salvaterra et al. 2008 considered two alter- 
native GC formation rates and found that these produce very 
similar results to those obtained by assuming that GC forma- 
tion follows the SFR. 
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(Co is a normalization constant). The best fit LF was de- 
rived by using the SHB rate -Rshb(z) for the two forma- 
tion scenarios (Eq. 1). Using these "best fit" LFs these au- 
thors calculated the predicted redshift distribution in the 
two cases and compared it with present sample of SHB 
redshifts (see Sect 2.2). By using the same procedure, we 
rederived the best fit parameters for the LF and deter- 
mined for the first time the local event rate, Rq (see Table 
1), independently in the two scenarios. We find that, by 
virtue of their longer r, dynamically formed systems would 
have to merge at a ^ 3 times higher rate than primordial 
systems in the local universe. 



The results in Table 1 are virtually insensitive to an in- 
crease of A2, whereas Lmin = L* / Ai cannot be increased 
by more than a factor of 10, without excluding the least lu- 
minous SHB detected so far (GRB050509B, Gehrels et al. 
2005). By decreasing Lmin = L*/Ai, the local SHB rate 
increases approximately as Rq cx (Lmin)" (Nakar 2006). 
Objects with luminosity Lmin observed by CGRO/BATSE 
or Swift /BAT at a limiting flux of -Flim, are detectable 
up to a maximum redshift of ^max which decreases with 
decreasing Lmin (for Swift/BAT F^m is ^ lphcm~^s~^; 
Sakamoto et al. 2008). Therefore SHBs with very low lu- 
minosity are above the limiting sensitivity of BATSE and 
BAT only in a fairly small volume, where they go off rarely. 
Current data do not constrain this part of the LF and the 
number of very weak SHBs (if any) remains unknown. 
It cannot be excluded that the volume afforded by the 
limiting sensitivity of BATSE and BAT is close to that 
required in order to sample the low-luminosity end of the 
LF of SHBs over the satellites' lifetime. However this ap- 
pears contrived and the low- luminosity end of the SHB LF 
might well extend below the BATSE and BAT sensitivity. 
We note in passing that the presently-known LF of long 
GRBs extends over more than 4 decades (see e.g. Liang et 
al. 2007, Guetta & Delia Valle 2007). In Sect 3 we discuss 
further the implications of a SBH LF that extends below 
the currently-estimated value. 



There exists by now evidence that SHBs are beamed 
in a relatively small solid angle. Similar to the case of long 
GRBs, the jet opening angle is inferred from the steepen- 
ing of the optical afterglow light curve. Fox et al. (2005) 
found a beaming factor of /^^ 50 for GRB 050709 and 
GRB 050724 (with /{, the fraction of the 47r solid angle 
within which the GRB is emitted). Soderberg et al (2006) 
inferred f^^ - 130 for GRB 051221A. We adopt a fidu- 
cial value of fi^^ ^ 100 and derive a beaming-corrected 
rate of po = f^^Ro ^ 130(/,-Vl00) and ^ 400(7^-7100) 
Gpc-'^yr-^ for SHBs originating in primordial and dy- 
namically formed binaries, respectively. The former esti- 
mate compares well with the lower end of the range es- 
timated by Belczynski et al. (2002) and Kalogera et al. 
(2004) for primordial DNSs ( 200 - 2800Gpc-3yr-i). 



2.2. The observed z distribution of SHBs 

In order to infer the redshift distribution of SHBs, we 
use the entire sample of presently known SHBs with red- 
shift determination, i.e. the sample of events in Table 2 of 
Salvaterra et al. (2008) plus GRB 071227, GRB 070429B 
and GRB 070714 at z = 0.384, z=0.902 and z=0.922, 
respectively (D'Avanzo et al. 2007; Berger et al. 2007b; 
D'Avanzo et al. 2008; Cenko et al. 2008). We excluded two 
of the bursts (GRB 061006, GRB 061201) from the sam- 
ple of Salvaterra et al. (2008) because uncertainties in the 
redshift determination □ We have verified that re- inserting 
these bursts in our sample does not affect significantly any 
of our results. 

In Figure 1 we compare the cumulative distribution 
of the observed redshifts with the expected cumulative 

distribution given by 

= ^f^^ fj^^ $o(L)dlogL , (3) 

where Lmax = A2L* = lOOL* and L,„ = 
max(LF, L*/Ai), with Lp the luminosity at z correspond- 
ing to Him- This is done in turn for the primordial and 
dynamical formation scenarios, by using the parameters 
given in Table 1 and Ai = 100 (see Fig. 1). A Kolmogorov- 
Smirnov (KS) test gives a probability of 0.04 and 0.2 
that the observed distribution is drawn from the distribu- 
tion expected for SHBs from primordial and dynamically 
formed NS-NS/BH binaries, respectively. These results are 
nearly insensitive to increasing values of Ai (see also in 
Sect. 2.1). 

We then considered the cumulative z-distribution that 
results from a combination of primordial and dynami- 
cally, formed NS-NS/BH systems. By varying their rel- 
ative incidence the closest matching with the observed 
z-distribution of SHBs was found for a ^ 60% contribu- 
tion from dynamically formed NS-NS/BH binaries (see the 
dot-dashed line in Figure 1). In order to get a rough lower 
limit on the incidence of dynamically formed NS-NS/BH 
systems we decreased their contribution until a KS prob- 
ability of 0.1 was reached: this gives a ~ 10% incidence of 
such systems. We note that the corresponding local SHB 
rate is i?o ^ 2.9 and 1.6 Gpc-^yr"! for a 60% and 10% 
fraction of dynamically formed mergers, respectively. 

Of course this incidence can be derived only approxi- 
mately, in consideration of the small sample of SHBs with 
redshift determination known so far. We emphasise that 
the relative contribution of the two populations of merg- 
ing binaries is determined primarily by the observed red- 
shift distribution of SHBs which can be constrained at 
present only through the 11 SHBs with known redshift. 

The above estimate compares well with the fraction 
(10 — 30%) of SHBs produced in GCs, as derived by 



^ However we note that D'Avanzo et al. (2008) have recently 
confirmed the redshift of GRB 061201. 

^ the other uncertainties in the expected cumulative distri- 
bution, see Eq. (3), play a negligible role here 
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Grindlay et al. (2006). We conclude that the present sam- 
ple of SHBs with redshift determination favors a bimodal 
origin of SHBs, with high-z bursts resulting primarily from 
the merging of primordial NS-NS/BH systems and low-z 
bursts produced mainly by dynamically- formed systems, 
as first proposed by Salvaterra et al. (2008). The contri- 
bution from the latter systems is unlikely to be < 10%. 
This conclusion might be altered if there existed a large 
population of high redshift (and thus mainly primordial) 
SHBs that has so far remained unidentified (Berger et al. 
2007a, O'Shaughnessy et al. 2008). 

On the other hand the association of a SHB with a 
GC would confirm unambiguously the dynamical forma- 
tion scenario. The large offset of some SHBs from their 
host galaxy is consistent with that of GCfl but direct ev- 
idence for this is still lacking. This issue may also be ad- 
dressed statistically by considering that, even though GCs 
are present in galaxies of virtually all typesH, most of them 
are in early type galaxies. Therefore SHBs from dynami- 
cally formed binaries, besides having on average lower red- 
shifts than those from primordial binaries, should be asso- 
ciated preferentially to early type galaxies. Since our sam- 
ple largely overlaps with that of Salvaterra et al. (2008), 
we obtain similar results to theirs in this respect, with 3 
in 3 SHBs associated to early type galaxies at z < 0.3, 
and 1 in 8 SBHs in association with early type galaxies at 
z > 0.3. While these results are in broad agreement with a 
bimodal origin of NS-NS/BH systems giving rise to SHBs, 
the low numbers involved prevents us from reaching any 
firm conclusion. 

The beaming-corrected SHB rate scales as 

Po-100(/,-Vl00)(Ai/100)"i?o (4) 

If dynamically formed DNS mergers represent 10% of 
the total SHB local rated, (Z^" VlOO)(^i/100)" ~ 0.1 
would be required for Eq. 4 to match the approximate 
merger inferred by Phinney et al. (1991) for GC DNSs 
{pdns ~ 3 Gpc~'^yr~^). This indicates that the rough 
estimate by Phinney et al. (1991) is about an order of 
magnitude lower than the beaming-corrected rate derived 
from SHB observations, a difference that we do not regard 
as crucial at this stage, in consideration of all uncertain- 
ties. 

3. Prospects for Gravitational Wave Detection 

Within the NS-NS/BH binary interpretation, the local 
rate of SHBs has clear implications for the number of 
merging events that can be detected with present and fu- 
ture ground-based GW interferometers. Based on the re- 
sults in Sect. 2, we estimate here the number of detectable 

* see Gehrels et al. (2005) for the case of GRB050509b and 
Bloom et al. (2006) for GRB060502b. 
^ except perhaps dwarf galaxies with Mv > — 13 
® Note that the contribution from BH-NS mergers is not in- 
cluded. 



GW events expected in Advanced LIGO/Virgo class in- 
terferometers due to SHBs originating from NS-NS/BH 

mergers, 

NgW ~ Wo[(l -gB)VDNS +gBVBHNs)] (5) 

' viooy V4Gpc-V~v Viooy 

^ I ^BHNS NT -1 
[(1 -.9b) +.gB-f7 )]yr , 

VdNSs 

where is the incidence of BH-NS systems among NS- 
NS/BH progenitors of SBHs, T^dns and Vbhns are the vol- 
umes corresponding to the horizon of GW interferometers 
for DNSs and BH-NS mergers respectively; the ratio be- 
tween these volumes is Vbhns/14)ns ^ 10. Here 77 = 1 for 
Advanced LIGO/Virgo and r/ = 3 x 10"* for LIGO/Virgo 
take into account the different volumes sampled by the 
two instrument classes (see e.g. Cutler & Thorne 2002). 
While it is estimated that only ~ 0.01 of primordial 
NS-NS/BH systems are BH-NS binaries (Belczynski et 
al. 2007), the fraction of such systems in GCs, though 
presently unknown, might well be higher, their main for- 
mation channel being the exchange interaction of a binary 
containing a NS and a cluster star with an isolated BH 
(Devecchi et al 2007). Dynamical star cluster simulations 
point to a fairly large population of BH-BH binaries in 
GCs, but have not yet provided estimates of (O'Leary 
et al. 2006). For instance for the globular cluster param- 
eters adopted by Sadowski et al. (2008) NS-NS/BH bina- 
ries are expected to form at a very slow pace, and it is not 
surprising that no such systems were produced in those 
simulations. 

We note that as long as gB^O.15 the rate of detectable 
GW events from BH-NS binaries will be higher than that 
from DNSs (see Eq. 5). We assume gB — 0.5 (and the 
other fiducial values in Eq.(5)) in the estimates of Nqw 
given below, and report also in parentheses the values cor- 
responding to gB ^ and 1, respectively. For Advanced 
LIGO/Virgo we find that events from primordial and dy- 
namically formed NS-NS/BH binaries are expected at a 
rate of 14 yr~^ and ~ 248 yr~^ (45, 450), respectively. 
For f^B = 0.5 the latter event rate is dominated by merging 
BH-NS binaries formed in GCs (~ 226 yr~^) and is close 
to the upper end of the range estimated by Kalogera et 
al. (2006) for primordial DNSs. For present generation in- 
terferometers Eq. 5 gives 1 event in 238 years and 1 event 
in 13 years (74, 7 years) from primordial and dynamically 
formed systems, respectively. 

There are presently substantial uncertainties in the 
values of the parameters in Eq.(5) and therefore a pre- 
cise estimate of the expected number of detectable GW 
event cannot be made yet. However it should be possible 
in the near future to constrain more tightly the uncertain 
parameters in Eq (5). Concerning the beaming factor, a 
few additional detailed studies of the optical afterglow of 
SHBs with currently available instrumentation will allow 
a more precise determination of f^^. Taken at face value, 
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present estimates range over a factor of ~ 2.6. The local 
SHB event rate Rq, while constrained by the large sample 
of SHBs revealed by BATSE , can vary up to a factor of 
~ 1.8 depending on the redshift distribution and incidence 
of dynamically formed SHBs (see Sect. 2.1 and 2.2). The 
latter is presently determined through 11 SHBs, but the 
sample of SHBs with secure redshift is steadily increasing 
and should triple by the end of the Swift mission. A more 
accurate determination of the lower end of the SHB lumi- 
nosity function will likely require more sensitive GRB de- 
tectors than currently available; however, as we discussed 
in Sect. 2.1 the lower limit Lmin will likely decrease (i.e. 
Ai will increase) resulting in a higher expected rate of 
detectable GW events. The incidence of NS-BH systems 
among dynamically formed SHB progenitors ge can be 
determined through more extensive dynamical GC simu- 
lations. The whole range of allowed values of gs {i-e. 0-j-l) 
translates into a factor of ~ 10 uncertainty in the number 
of GW events in Eq.(5). 

Our treatment assumes that SHBs with similar prop- 
erties are produced both by DNS and BH-NS mergers. 
From the point of view of the models, this appears to be 
a reasonable ansatz as the torus of matter at nuclear den- 
sities is expected to originate from the lightest collapsed 
object in the system and provide the main source of en- 
ergy for the SHB. Moreover, observations have not yet 
shown clear evidence that SHBs comprise different sub- 
classes. Alternatively, if all SHBs detected so far came 
from DNSs, the total GW event rate from dynamically 
formed DNS and BH-NS mergers would be higher (unless 
qb = ^ holds). On the more pessimistic side, if a fraction 
h of SHBs were due to an entirely different phenomenon 
[e.g. giant flares from Soft Gamma Repeaters, Hurley et 
al. 2005) all the rates given above would be decreased by 
a factor ~ (1 — h). 

According to Sadowski et al. (2008), the rate of BH-BH 
mergers formed in globular clusters is high and can give 
rise to 0.01-1 event per year in LIGO and 25-3000 events 
per year in Advanced LIGO, a major uncertainty being 
the initial stellar mass fraction in dense clusters. Other au- 
thors envisage substantially lower rates (Portegies Zwart 
and McMillan 2000, O'Leary et al. 2006). Once the ratio 
between NS-NS /BH and BH-BH binaries in globular clus- 
ters will be better known from simulations, the merging 
rate of dynamically formed NS-NS/BH binaries inferred 
from SHBs can provide an independent normalization for 
the BH-BH merger rate. 

We note that an interferometer's horizon would in- 
crease by a factor of ~ 2.4 for a GW detection in co- 
incidence with a SHB detection (Cutler & Thorne 2002). 
This is because of the higher sensitivity afforded by having 
an independent knowledge of the time of occurrence and 
sky position of the event. The above factor takes also into 
account that, if SHBs are beamed along the binary's an- 
gular momentum axis, then the SHBs that are detected 
at the earth are also the ones whose orientation angle 
maximises the GW signal. The larger horizon translates 
into a factor of ~ 15 increase in both Vdns and Vbhns- 



Considering that only a fraction ~ 0.01 of the SHBs 
are beamed toward us, the coincident event rate would be 
~ 0.15(/^"Vl00) of the rates derived from Eq. 5. 

4. Conclusions 

In this paper we have presented a simple method for in- 
ferring the merging rate and corresponding rate of de- 
tectable GW signals from coalescing NS-NS/BH binaries. 
The method is based on the assumption that at least a 
sizable fraction of SHBs originates in the merging of such 
binaries and exploits the most recent results from observa- 
tions to infer the luminosity function, redshift distribution 
and degree of beaming of SHBs. We find further evidence 
in favor of a bimodal origin of SHB progenitors, with the 
merging of primordial binaries dominating at high red- 
shift, whereas at lower redshift at least ^ 10% (and prob- 
ably about half) of the events arises from coalescing bina- 
ries formed dynamically in GCs. We have shown that the 
latter events make the expected local merging rate higher. 

The accuracy of the GW event rate that can be esti- 
mated at present is hampered by our present knowledge of 
the parameters in Eq.(5). Maximmn uncertainties in the 
SHB beaming factor (/b) and local rate (i?o) add up to 
a factor of ~ 5, whereas the ignorance of the incidence 
qb of BH-NS binaries among dynamically formed systems 
adds a factor of 10 uncertainty. This results in total uncer- 
tainty by a factor of ~ 50 in the detectable GW event rate 
in Eq. (5). Though large, this is not distant from the un- 
certainties of other estimates of GW events in literature. 
For instance Kalogera et al. (2004) estimate the galactic 
double neutron star merger rate with an uncertainty of 
factor of 14. Similarly, the results of the dynamical sim- 
ulations of Sadowski et al. (2008) for BH-BH mergers are 
affected by a 1.5 decade wide uncertainty in the initial 
mass fraction of globular clusters. We expect that further 
SHBs observations in Swift era will lead to a more ac- 
curate determination of fb and Rq, while more advanced 
dynamical globular cluster simulations will decrease the 
range of allowed values of gs. Progress in this area might 
well take place at a sufficiently fast pace, that our method 
can yield accurate predictions before GW events from co- 
alescing binaries are detected in large numbers. 

We have shown that a sizeable fraction of detectable 
GW events is expected to be coincident with SHBs (de- 
pending on the SHB beaming factor): this provides a new 
interesting perspective for the Advanced LIGO /Virgo era. 
We conclude that SHBs can provide a powerful means of 
inferring the GW event rate from coalescing binaries con- 
taining two collapsed objects. 
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Table 1. Model parameters and local rates 



DNSs (BH-NS) 


L* 

[10^^ ergs^^] 


a 




Ro Rate(z=0) 
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Fig. 1. Observed (histogram) and expected (curves) cu- 
mulative 2-distribution functions for the different time- 
delay probabilities expected in the two scenarios. Model 
distributions are shown for different combination of pri- 
mordial and dynamically formed NS-NS/BH systems. 
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